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ABSTRACT 



In this work we analyze the physical properties of a sample of 56 spectroscopically selected star-forming (SF) Lyo- emitting galaxies 
at 2.0<z<3.5 using both a spectral energy distribution (SED) fitting procedure from rest-frame UV to mid-IR and direct 160^m 
observations taken with the Photodetector Array Camera & Spectrometer (PACS) instrument onboard Herschel Space Observatory. 
We define LAEs as those Lyo- emitting galaxies whose rest-frame Lycr equivalent widths (Lyo- EW resl _f rw „ e ) are above 20A, the typical 
threshold in narrow-band searches. Lyo- emitting galaxies with Lyo- EW M _/ rame are called non-LAEs. As a result of an individual 
SED fitting for each object, we find that the studied sample of LAEs contains galaxies with ages mostly below lOOMyr and a wide 
variety of dust attenuations, SFRs, and stellar masses. The heterogeneity in the physical properties is also seen in the morphology, 
ranging from bulge-like galaxies to highly clumpy systems. In this way, we find that LAEs at 2.0<z<3.5 are very diverse, and do 
not have a bimodal nature, as suggested in previous works. Furthermore, the main difference between LAEs and non-LAEs is their 
dust attenuation, because LAEs are not as dusty as non-LAEs. On the FIR side, four galaxies of the sample (two LAEs and two 
non-LAEs) have PACS-FIR counterparts. Their total IR luminosity place all of them in the ULIRG regime and are all dusty objects, 
with A 12 oo^4mag. This is an indication from direct FIR measurements that dust and Lyo- emission are not mutually exclusive. This 
population of red and dusty LAEs is not seen at z~0.3, suggesting an evolution with redshift of the IR nature of galaxies selected via 
their Lya emission. 

Key words, cosmology: observations - galaxies: stellar populations, morphology, infrarred. 



1. Introduction 

Over the past years, a large number of galaxies have been found 
at different redshifts, from the local universe up to z~7 and 
even beyond. One of the most successful selection methods for 
searching for high-redshift galaxies is the narrow-band tech- 
nique. It employs a combination of narrow and broad band filters 
to isolate the Lya emission in the spectrum of a galaxy and to 
constrain its nearby continuum, respectively. This method segre- 
gates galaxies with a Lya emission whose rest-frame Lyo- equiv- 
alent width (Lya EW rejr _/ rame ) is typically above 20A. They are 
called Lye emitters (LAEs). 
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* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



Many efforts have b een aimed at looking for LAEs in a 
wide r ange o f redshifts (|Deh arveng et al., 200 8c ICowie 
20101 I201 U iBongiovanni et all l2010b iGuaita et al 
Hul 
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of LAEs have been mostly analyzed by fit ting their observed 
spectr al energy distribution (SED) with iBruzual & Chariot! 
(2003, hereafter BC03) templates. This allows the determination 
of their d ust attenuation, star-formation r at e (SFR), age or stel 
lar m a ss jFinkelstein etail 120081 [2009alfbl: iNiisson et all 12007 
2009, 20 1 It ICowie etall 120 lit lOno etall 12010c iLaieTal 
20081 iPirzkal et all 12007c IGawiser et all 1200 



2007). 



BC03 templates do not take the dust emission in the FIR 
into account and, even if they did, the lack of FIR information 
for LAEs, mainly at z>2, would produce important physical pa- 
rameters, such as dust attenuation and SFR, suffering from large 
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uncertainties. In lOteo et all (1201 lalbl) we look for mid-IR/FIR 
counterparts of a sample of LAEs at z~0.3 by using PACS- 
lOOyum, PACS-160//m, and MIPS-24/mi data, finding that a high 
percentage of them (~75%) are detected at those wavelengths. 
These detections enable us to determine their IR nature, dust at- 
tenuation, and SFR without the intrinsic uncertainties of SED 
fitting. As a result, we find that LAEs at z~0.3 are among the 
least dusty galaxies at that redshift and the majority have total 
IR luminosities below lO n L . 

The 2<z<3.5 range (1.7 to 3.2 Gyr after the Big Bang) 
is where th e SFR density of the Universe is at a maximum 
dHogg et all 119981: lHopkinsl l2004t iHopkins & Beacoml [2006; 
iPerez Gonzalez et all 120051 iLe Floc'h et all 120051) . and, with 
respect to LAEs, i t has only started to be studied recently. 
INilsson et alJ (2009) demonstrate that there is a significant evo- 
lution in the physical properties of LAEs between z~3.0 and 
z~2.3, with a spread in their SEDs which is greater at z~2.3 than 
at z~3.0. This indicates that LAEs at z~2.3 are more massive, 
older, and/or dustier than those at higher redshifts. INilsson et alJ 
(1201 lh find that both the dust attenuation and stellar mass of 
LAEs at z~2.3 are high, mainly spanning Ay- 0.0-2.5 mag and 
log(M»/M )= 8.5-1 1.0, respectively. They also find that physical 
proper ties of LAEs at that redshift are very diverse. Guait a et"afl 
(2011) obtain robust determinations of mass and dust attenua- 
tion, log(M»/M ) = 8.6[8.4-9.1] and E(B-V) = 0.22[0.00-0.31], 
for LAEs at z~2.1. Furthermore, comparing with a sample of 
LAEs at z~3.1, they find that LAEs at z~2.1 tend to be dustier 
and show higher instantaneo us SFR than those at z~3.1, and 
the properties are also diverse. iBond et all d201 ll) claim between 
z~3.1 and ~2.1 for an evolution of the morphological proper- 
ties of LAEs, with the median half-light radii rising with de- 
creasing redshift. They also report that LAEs at z~2.1 are bigger 
for galaxies with higher stellar mass, star formation rate, and 
dust obscuration. Therefore, 2.0<z<3.0 is an interesting redshift 
range for several reasons: 1) large samples of LAEs are being 
collected, 2) in this redshift range, LAEs tend to have bright ob- 
servable fluxes, which allows a multiwavelength coverage with 
a higher signal-to-noise ratio than at higher redshifts, 3) it is a 
redshift range where the physical properties of LAEs seem to be 
changing significantly and become very diverse, while this has 
not been clearly reported at higher redshifts. 

In this work, we focus on a sample of spectroscopi- 
cally selected Lya emitting galaxies at 2<z<3.5. The spec- 
troscopic selection segregates galaxies with Lya EW rest-frame 
either above or below 20A, the typical threshold in narrow- 
band searches. To place our contribution in a common base 
with narrow-band selected galaxies, we distinguish the follow- 
ing throughout the work: LAEs: Lya emitting galaxies with Lya 
EW J - fJ ,_/ rame ,>20 A; non-LAEs: Lya emitting galaxies with Lya 
EW J - ej ,_/ rame ,<20 A. The main objective of this work is to give 
a first glimpse at the FIR properties of LAEs at 2.0<z<3.5 by 
using deep FIR data coming from H e rsche l-PACS observations 
(Pogl itschetafl 120101: iPilbratt et all 120101) . Furthermore, with 
the aim of better understanding this population in that redshift 
range, we also study their physical properties by employing a 
SED fitting procedure with BC03 templates. This allows the 
comparison between SED fitting and IR-based results. 

This paper is divided in two main parts. First, in Sections [2] 
[3] |U and [5] we describe the optical to mid-IR data used in this 
work, perform SED fittings with BC03 templates, and carry out 
a morphological analysis. Then, in Section [6] we focus on the 
FIR side of their SED, studying total IR luminosities, IR nature, 
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Fig. 1. Redshift distribution of the sources in the final sample. 
We study objects at 2 < z < 3.5. 

total SFR, and dust attenuation for the FIR-detected galaxies. 
Section|7]gives the conclusions of this work. 

Throughout this paper we assume a flat universe with 
(Cl m , £2a, ho ) _ (0.3, 0.7, 0.7). A ll magnitudes are listed in the 
AB system dOke & GunnUl983l) . 

2. Optical data and object selection 

The sample used in this work was selected from the final data 
release of t he VIMOS spectroscopic campaign in the GOODS- 
South field dPopesso et all l2009t iBafestra et alll2010l) . VIMOS 
spectra have a quality flag assigned as follows: A (secure clas- 
sification), B (likely classification), and C (tentative classifica- 
tion). From all the spectra in the survey, we only selected those 
galaxies that exhibit a Lya emission in their spectrum and which 
are flagged as A or B. This yields a sample of 144 objects. The 
wavelength coverage of the spectra implies that the LAEs are 
distributed with in 2.0<z<3.5. Thro ughout this work, the objects 
are named as in lPopesso et alJ d2009l) . 

With the aim of carrying out SED fittings for the studied 
galaxies, we used the MUSIC multiwavelength photometric cat- 
alog (from 0.3 to 24.0 /urn) of a large and dee p area in the 
GOO DS-South Field. Its first version was made by Grazian et al. 
(2006), and a subsequent update was elaborated bv lSantini et al. 
(2009). The former uses F435W, F606W, F775W, and F850LP 
ACS images, JHKs ISAAC data, mid-IR data provided by IRAC 
instrument (3.6, 4.5, 5.8, and 8.0 //m), and publicly available U- 
band data from WFI and VIMOS. The z-band of ACS GOODS 
frames and Ks of VLT images were used to select objects in the 
field, yielding a unique, autoconsistent catalog. The details of 
the catalogs (object detection s, PSF-matchings, l imiting magni- 
tudes, etc.) can be found in iGrazian et alJ d2006l) . The second 
version includes objects selected from the IRAC-4.5 fim image 
and, therefore, sources detected in that wavelength but very faint 
or undetected in K s band. MIPS-24yum photometry was also in- 
cluded, although it will not be considered in the SED fittings 
(see Section [3]). Owing to these improvements, we adopted the 
second version of the photometric catalog. 
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Fig. 2. Histogram of the magnitudes in the V band of the final 
sample of 56 LAEs. The values appearing here are similar to 
those in the UV-bright sample of Guait a et ail (1201 lh . 

To build our final sample of Lya emitting galaxies, we 
matched the coordinates of 144 spectra with the MUSIC cata- 
log, resulting in a total of 70 objects. The loss of 50% of the 
sources is explained by the area covered by the VIMOS spec- 
troscopic observations being larger than the GOODS-MUSIC 
footprint. From that common sample, we visually inspected each 
spectrum, ruling out objects with no clear Lya emission or with 
AGN ionization lines in their spectra. We finally ended with a 
clean sample of 56 Lya emitting galaxies whose optical spec- 
tra suggest that they have a star-forming (SF) nature. Figure [1] 
shows the redshift distribution of the objects in the final sample. 
It can be seen that most of them are at 2.0<z<2.7, although we 
also have a significant number of sources at z>3. Figure|2]shows 
the distribution of the observed V-band magnitude of the stud- 
ied galaxies. Most of them are brighter than 25.5 mag in that 
band, which is comp arable to the bright ness of the UV-bright 
subsample defined in iGuaita et al.l d201 lh . The main advantage 
of working with a continuum-bright sample is that it is possible 
to carry out individual S EP fitting for each ga laxy, avoiding the 
uncertainties of stacking (iNilsson et alll201 ll) . 

3. SED fitting method 

We performed SED fittings with BC03 templates for the sam- 
ple of 56 galaxies by using th e Zurich Extraga l actic Bayesian 
Redshift Analyzer (ZEBRA, iFeldmann etaU l2006h . In its 
maximum-likelihood mode, ZEBRA employs a \ 2 minimiza- 
tion algorithm over the templates to find the one that fits the 
observed SED of each object best. In this process, we excluded 
the MIPS-24//m measurements, since the fluxes in this band 
have a significant contribution of warm dust and PAH molecules 
emission, which are not taken into account in the BC03 tem- 
plates. We employed GALAXEV, which is provided by BC03, to 
build a large sa mple of templates. In this process we adopted a 
Salpeter ( 1955) initial mass function, distributing stars from 0.1 
to 1OOM . This IMF has mo stly been used in previous works 
in the studied redshift range dGuaita et all 1201 it INilsson et all 



IL~ai et a 120081: lOno et all 120101: iGawiser et all [2006, 
2007). Since metallicities do not significantly modify the shape 
of the SEDs of galaxies, we adopted Z=0.4 Z o . Different au- 
thors employ different values: lLai et al.l d2008h and IGuaita et al.l 
d201 lh consider the solar metallicity, whereas [Ono et al.T T2010) 
utilize a value of O.2Z . To check the validity of our choice, 
we ran ZEBRA with BC03 templates associated to different 
metallicities and we find similar results for the other parameters, 
within the uncertainties. This, at the same time, is an indication 
that metallicity is difficult to determine accurately with SED fit- 
ting. For the SFR, we use mode ls with temporally constant val- 
ues, as in many previous works ([Gawiser et ail 120061: lLai et all 
120081: lOno etalll2010t IGuaita et all 1201 lb iKornei et all l2010h . 
In this case, the SFR can be obtained from the rest-frame UV 
fluxes once the templates are norma l ized t o the observed pho- 
tometry and employing the K ennicuttl (fl998 ) calibration between 
SFR and the UV luminosity. We also tried to use exponentially 
decreasing SFR, but the differences betw een the results fo r age 
and dust attenuation were not significant. lOno et~ail feOlO) also 
compare the results obtained with constant and varying SFR and 
find similar results in both cases, which indicates that the history 
of the star formation in a galaxy is hard to constrain with SED- 
fitting methods. Furthermore, adopting an exponentially varying 
SFR would add a new degree of freedom in the process, the star- 
formation time scale, introducing more uncertainties in the de- 
termination of the other parameters. 

The ages considered here span from 1 Myr to 3.4 Gyr, the age 
of the universe at the minimum redsh ift of the sample, z ~2.0. 
Dust attenuation is included via the Calzetti et a 1. (2000) law 
with values of the color excess in the stellar continuum, E(B- 
V) s , ranging from to 0.7 in steps of 0.05. We also inclu de inter- 
galact ic medium absorption adopting the prescription of Madau 
d 19951) . Stellar masses are obtained from SFR and age, accord- 
ing to the assumed temporal variation of the SFR. Regarding the 
uncertainties in the fitted values, ZEBRA provides a parameter 
related to the probability that one template truly represents the 
observed photometry of a given source. In the SED fitting of 
each galaxy, we select all the templates whose associated proba- 
bility is above 68% and consider the uncertainty of each param- 
eter as the average of the values associated to that parameter on 
the selected templates with more than 60% of probability. 



4. Stellar populations 

Figure [3] and Table I A. 1 1 show the results of the SED fitting for 
the studied galaxies. The dust attenuation of the studied LAEs, 
parametrized by the color excess in the stellar continuum E(B- 
V) s , ranges from to 0.3, with a median value of 0.15. The 
uncertainties in dust attenuation have a lower limit equal to 
the sampling of this parameter in the templates, that is, AE(B- 
V)j~0.05, and are also affected by the well-known degener- 
acy between age, dust attenuation, and metallicity of the SED- 
fitting methods. Despite the uncertainties, the values obtained 
here are distributed within approximately the same interval as in 
IGuaita et al.l d201 lh for the UV-bright sample, although their me- 
dian value, E(B-V)~0.32, is higher than that found in the present 
work. 

Regarding age, the studied LAEs are young galaxies, most 
of them having ages below 100 Myr. However, there are some 
LAEs, ~20%, whose SEDs are compatible with older popula- 
tions, above 300 Myr. Although the threshold age in the gener- 
ation of the templates was set to the age of the universe at z~2, 
there is no violation of the age of the universe for galaxies at 



3 



I. Oteo et al.: Stellar populations and Herschel FIR-counterparts of Lyo- emitters at 2 < z < 3.5 in GOODS-South 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 
E(B-V) 



500 1000 
Age [Myr] 



1500 




20 40 60 80 100 120 
SFR UV [M yr" 1 ] 



Log(M,/M 3 ) 



Fig. 3. Distribution of the physical properties for the sample of 56 LAEs. Blue and red-shaded histograms represent LAEs and 
non-LAEs, respectively. Inset plots are more detailed representations of the zones of the histograms where most objects are located. 



higher redshifts. The uncertainties in age are typically below the 
20%, which hints towards robust age estimations. 

Intense rest-frame optical emission lines, such [Oil], [OIII] 
or Ho-, could have a strong effec t on the observed broadband 
photometry of a galaxy. At z~0.3. ICowie et all d201 lb employed 
their rest-frame optical spectroscopic observations of their LAEs 
to analyze the influence of rest-frame optical emission lines on 
the derived parameters in a SED fitting procedure. As a result, 
it was found that not subtracting those lines overestimates of the 
ages of the LAEs in their sample. In our work, given that we 
do not have spectral information about the [Oil], [OIII], or Ho- 
emission for the studied galaxies, we can not correct for this 
effect, and, so the ages obtained here should be considered as 
an upper limit of the real values. 

Since the studied LAEs have a bright UV continuum, 
their dust-uncorrected SFRs, SFR(/y„„ eorrf , c/ ^, are high, rang- 
ing from about 2 to 100 M yr _1 , with a median value of 
S¥Ru VtUncorrecte d~\0 M yr _1 . These values are slightl y higher 
than those reported at lower redshifts. In lOteo et alJ (1201 lbl) 



we studied the SFRuv,uncorrected f° r LAEs at z~0.3, finding that 
most of them have SFRuv,uncorrected below 5 M yr _1 . In the 
present work we do not see such low values because of the UV- 
brightness of the studied LAEs. However, the high values found 
at 2.0<z<3.5 are not seen at z~0.3. This could suggest that the 
upper limit of the SFRuv tmc0 rrecte4 m LAEs is increasing from 
0.3 to z>2, or similarly, that there is an evolution in the upper 
limit o f the UV luminos ity, Lyy. The surveyed comoving vol- 
ume in lOteo et all (1201 lbl) and in the present work are similar 
given that the GALEX observations used in the former cover a 
much larger area of the sky than those used in the present work 
because the large r surveyed area at z~ 0.3 makes up for the differ- 
ence in redshifts. ICowie et al.ld201 lb report a strong evolution in 
the Lya luminosity function of LAEs between z~0.3 and ~1.0, 
which could be related to the evolution in the upper limit of Lj/y 
reported above. This suggests that there is a noticeable change 
in the UV properties of galaxies selected via their Lya emission 
from z-0.3 to z>2.0. 



4 



I. Oteo et al.: Stellar populations and Herschel FIR-counterparts of hya emitters at 2 < z < 3.5 in GOODS-South 




-0.2 0.0 0.2 0.4 0.6 0.8 

V-I [mag] 

Fig. 4. Color-color diagram containing the locus of LBGs according to lPentericci et al.l (l2010h . Filled and open dots represent LAEs 
and non-LAEs, respectively. Orange and blue dots represent galaxies at z>2.75 and z<2.75, respectively. Red dots are PACS- 
detected galaxies and shaded zone is the area where LBGs are located. For clarity, the galaxies with a U-V color greater than 5.5 
are assigned U-V=5.5. 



The stellar mass of the studied LAEs spans within 
7.O<Log(M»/M )<9.5. More than 70% of the studied Lya emit- 
ting galaxies are detected i n IRAC-3.6^m u nder a similar limit- 
ing luminosity to the one in lLai et aT . (2008). In this way, we are 
working with a sample of IRAC -bright sources, and they should 
be among the most massive LAEs in the studied redshift range. 
The range o f stellar mass es obtained here is compatible with 
the results of Gua ita et al. (1201 ll) for the UV-bright sample and 
with t hose at z~3.0 of lLai et all (120081) and lGawiser et all (120061 
120071) for both IRAC-detected and IRAC-undetected stacked 
samples. 

After inspecting the individual SED fittings we conclude that 
there is a wide range of shapes, from almost-dust free cases, 
with a negative UV continuum slope, to highly attenuated LAEs, 
with a positive UV continuum slope. The rest-frame optical-to- 
UV colors also show a significant variation from blue to red ob- 
jects. This way, analyzing individual SEDs for each object, we 
find LAEs with very different properties, rather than the dou- 
ble nature suggested in previous works, which could be t he re- 
sult of the stacking metho d being e mployed ( Lai et ail [2008 : 
iGawiser et al 1 I2006L 120071) . Actually. iNUsson et al.l (1201 1 stud- 
ied the effects of stacking in their sample of LAEs at z~2.3, find- 
ing that, while stellar masses are robust to stacking, ages and 
dust attenuation tend to be incorrectly determined. Therefore, it 
is clear that to analyze the physical properties of LAEs at high 
redshift, individual fits to the observed photometry are needed, 
although it is challenging for those LAEs at z>3 owing to pho- 
to metric hmi t ations. Actu ally, only a few LAEs of the samples 
of lOno et al] (1201 Ol) and lLai et al.l (120081) are detected in the 



Ks band, which, at their redshifts, is quite important for sam- 
pling the Balmer break and obtaining accurate values of age 
and mass. Results based on stacking analysis must be taken with 
care, most importantly in those cases where the SEDs are known 
to have very diverse shapes and propertie s, as shown i n this 
and previous works to happen at 2.0<z<3.5 dGuaita et al.L 1201 it 
iNilsson etal.Ll2011l) . 

As mentioned in Section [1] the spectroscopic selec- 
tion used in this work allows isolating sources with Lyar 
EW re «_/ rflm e<20 A (non-LAEs). However, the number of such 
sources in our sample is low. Non-LAEs are also represented 
in Figure [3] It can be seen that the age, SFR and mass of 
these galaxies are in the same range as those for LAEs. The 
main difference is the dust attenuation: objects with an attenu- 
ated Lya emission are among the dustiest objects in the sam- 
ple. This points towards the idea that LAEs are among the 
least dusty galaxies at each redshift , as reported in other wor ks 
for different redshifts (ICowie et all 1201 lUOteo et al.L l2l)TTallbt 
iPentericci etafl 120071 120 1 Ot iKornei etaiT 120 1 Oh . 

Figure |4] represents the locus of Lyman Break Galaxies 
(LBGs) at z~3 .0 fo r the filters in the MUSIC survey 
(iPentericci et al.Ll2010l) . It can be seen that most LAEs at z>2.75 
could have been also selected as LBGs. Most LAEs at z<2.75 
do not meet the LBG criterion, because for those redshifts, the 
Lyman break is not located between the U and V bands, but in 
bluer wavelengths. The subsample of LAEs with colors of LBGs 
(LAE-LBGs) is the bridge population between these two kinds 
of galaxies. The principal difference between LAE-LBGs and 
the total sample of LAEs is the dust attenuation: LAE-LBGs 
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tend to be dustier than other LAEs. In Figure |H this tendency 
is reflected in their V-I color. At our redshifts, this color mea- 
sures the UV continuum slope, and higher values correspond to 
greater dust attenuation. LAE-LBGs have higher V-I color (red- 
der UV co ntinuum) than ot h er LA Es of the sample, since they 
are dustier. iPentericci et al.l d2010h find that LBGs with a Lya 
emission at our redshift have lower dust attenuation than those 
without a Lya- line. The combination of both results points to- 
ward a sequence of increasing dust attenuation, from the less- 
dusty non-LBG LAEs, to dustier LAE-LBGs, and finally those 
LBGs without Lya emission, so they could be the same popula- 
tion differentiated by dust attenuation. 



5. Morphology 

By using HST/ACS images in GOODS-South we carried out a 
morphological analysis of the 56 studi ed Lya emitting galaxies . 
We adopted the classification given in Elmegreen et al. (2009), 
and visually classify them into three groups: 37 clump clusters 
(CC), 8 chain galaxies (CH), and 1 1 spiral-bulges (SPB). This 
classification was done by three different people and the results 
were similar in the three cases. There is a wide variety of sizes 
and morphologies, from highly compact to highly clumpy ob- 
jects. This heterogeneity in morphology is compatible with the 
wide range of physical properties for our LAEs reported in the 
previous section, indicating that galaxies segregated by their Lya 
emission at our redshift do not have specific properties, but in- 
stead can be very diverse. 

We fit the SPB LAEs, whic h are all at z<2.75 , with lSersid 
(119681) profiles by using GALFIT dPeng et alll2010h in the z-band 
images. At our redshift, this band represents the emission of the 
rest-frame optical light near the Balmer break. The Sersic pro- 
files can be analytically described by 



I(r) = I b (0)exp 



-b»(r/R eS ) 



l/n 



(1) 



where h(0) is the central intensity, 7? e g- the effective radius, and 
n the Sersic index. Figure [5] shows how GALFIT performs the 
fittings in three randomly selected SPB LAEs, which are repre- 
sentative of the behavior of th e total sample of SPB LAEs. It can 
be seen that the lSersicl {l968) profiles describe the morphology 
of these objects quite well. To obtain the intrinsic R e ^ in kpc 
from the output value given by GALFIT (in pixels), we made use 
of the assumed cosmology and the ACS pixel scale. As the re- 
sult, we find that the SPB LAEs have effective radii ranging from 
2.6 to 3.1 kpc. The uncertainty in the effective radius, which is 
provided by GALFIT, is less than 5%. This range is compatible 
with the tail of th e distribution at large effective radii found in 
Bon d et alJ (l2009h at z~3.1. in lBondetafl (l201lh at z~2.1, and 



Oteo et al. (1201 lbl) at z~ 0.3. Therefore, an evolution in the 



physical sizes of LAEs with redshift is not seen with the present 
data. 

Regarding the Lya EW„„_/ ra/K , we do not find any sig- 
nificant difference in the morphology of LAEs and non-LAEs, 
because of t he scarcity of gala xies in the non-LAE sample. 
At z~0.3, in lOteo et all d2011bh we find that LAEs and non- 
LAEs have a clear difference in their morphology, both in shape 
and size, which indicates that at that redshift Lya photons are 
escaping from small a nd irregular/merging galaxies. However, 
Pentericci et alJ (1201 Oh show that at our redshift, the rest-frame 
UV morphology does not depend strongly on the presence of the 
Lya emission. 







n 















Fig. 5. Examples of GALFIT fittings for three SPB LAEs. 
Left: Images of the LAEs. Middle: Sersic fitted profiles. Right: 
Residuals of the fittings. 



6. Herschel FIR counterparts 

GOODS-South was observed with PACS-70yum, PACS-lOOyum, 
and PACS-160/im in the frame of the PACS Evolutionary Probe 
project (PEP, PI D. Lutz). PEP is the Herschel Guaranteed 
Time Key-Project designed to obtain the best profit from 
Herschel instrume ntation to study the FIR galaxy population 
dLutz et al.L 1201 lb . PACS fluxes used in this work were ex- 
tracted using MIPS-24/mi position priors with, at least, a 3<x 
significance. Limiting fluxes in PACS-70/mi, PACS-100/im, and 
PACS-160/im are l.OmJy, l.lmJy, and 2.0mJy, respectively. 
We look for possible FIR counterparts of our galaxies within 
2", which is the typical astrometric uncertainty in the position 
of the sources, finding four detections in at least the PACS- 
160//m band. These counterparts are direct evidence of dust 
emission in high-redshift SF sources. As seen in Figs. ??, ?? 
and ??, GOODS.LRb.001_q2_9.1, GOODS_LRb_001.q3_9.1, 
and GOODS_LRb_001_ql_8_l are isolated sources and, there- 
fore, the dust emission is coming entirely from these sources. 
However, GOODS_LRb_dec06_3_q3_60_l has a nearby galaxy 
whose emission in the FIR could contaminate the FIR flux of 
this galaxy. Two of the PACS-detected galaxies are within the 
LAE group and the other two have Lya EW rf .„_/ rame compatible 
with their being non-LAEs. Figure [6] represents the rest-frame 
UV to FIR SED of the four PACS-detected galaxies. 

Table Q] shows the dust attenuation and SFRuv, u „ C orrected 
of the four PACS-detected galaxies, as derived from SED 
fitting. Their rest-frame UV to mid-IR SEDs are compati- 
ble with high values of dust attenuation, and the two LAEs 
also exhibit greater SFRuv.imcorreaed than the median value 
of the sample. Their large SFRuv.uncorrected indicates that 
they have high value of UV luminosity, compatible with 
the greater li kelihood tha t UV-b right galaxies are detected 
in the FIR (iReddv et all |2010). GOODS_LRb_001.q2_9.1, 
GOODS_LRb.001_q3.9_l, and GOODS_LRb_dec06_3_q3.60_l 
are young objects with masses similar to the median population. 
However, GOODS_LRb_001.ql.8_l is an old and dusty LAE, 
and it is among the most massive objects of the sample. This 
is one of the most interesting galaxies in our sample, since its 
physical properties are opposite the classical idea of an LAE, 
i.e. a young, less massive, and dust-free galaxy. Three out of the 
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four PACS -detected galaxies would have been selected via the 
Lyman break technique (see Figure|4|l. The other PACS-detected 
source would have not been segregated as LBG either, since its 
redshift is not high enough to place the Lyman break between 
the U and V filters. It can be seen that PACS-detected sources 
have among the highest V-I colors of the total sample, indicating 
that they have a high-attenuated UV continuum, compatible with 
being among the dustiest galaxies in our sample. 

6.1. IR luminosities 

Accurate values of the total 8-1000pm IR luminosities, Ljr, can 
be obtained for the PACS-detected galaxies. Since we do not 
have full coverage of the dust emission peak, we cannot calcu- 
late Ljr by direct integration of th e FIR SEP. Ins t ead, w e con- 
vert PACS fluxes into L IR by usin dCharv & Elbazl (l200ll) (here- 
after CE01) templates, which have only one solution for L/r for 
each FIR flux and redshift. Obtaining L/« by using single FIR 
band extrapolations has been employed with suc cessful results 
in previous works in galaxies a t simila r redshifts (|Nordon et al.L 
2010; Elbaz etalll201(ll2011h . In fact lElbaz et al.1 (1201 Oh com- 
bine PACS and SPIRE measurements to show that fitting CE01 
templates to the PACS-160yum flux provides reliable estimates 
of L/R at our redshift. The typical uncertainties of the derived 
values are related to the accuracy of single band extrapolations 
rather than to the uncertainties of the PACS fluxes t hemselves, 
and are typically less than 0.2 dex (Elb az et al 1 1201 lh . Actually, 
for galaxies in our redshift range, PACS-160yi/m band extrapo- 
lation s give the most ac curate results for single band extrapola- 
tions (|Emizetifl[20i3). 



6.2. LAE-ULIRGs relation 

Previous results suggest that some LAEs at our redshift 
range are also ultraluminous infrared ga l axies (ULIRGs) 
(iNilsson & Moiled. l2009t iNilsson et all 1201 lb lOno et all l2010t 
Chapman et all 120051) . but direct detections in the FIR around 
the dust emission peak have not been reported yet. In the present 
work, according to their L/«, the four PACS-detected galax- 
ies have an ULIRG nature. Owing to the depth of the PACS 
observations used in this work, the limiting L/« is about 10 12 
L at z~2.2. Therefore, there could be some PACS-undetected 
LAEs in our sample, the dustiest at the highest redshifts, that fall 
in the ULIRG class. Actually, given that FIR detections segre- 
gate dusty galaxies, it is expected that the dustiest LAEs, E(B- 
V)>0.2 of the sample have an ULIRG nature even when PACS- 
undetected. Thus, more LAEs with an ULIRG nature should be 
expected. The number of LAEs with an ULIRG nature f ound, 
although scarce, is very important. In lOteo et all d201 1 alibi) we 
studied the FIR properties of a sample of 23 mid-IR/FIR detected 
LAEs at z~0.3, finding that most LAEs are in the normal SF 
galaxy regime, L/R<lO n L , with only one having L/r~10 1l5 L o , 
and none above that value. The discovery of LAEs with an 
ULIRG nature at z~2.4-2.8 suggests that the IR properties of 
LAEs evolve from z~0.3 to ~2.5 in the sense that there is a pop- 
ulation of red and dusty LAEs that is not seen at lower redshifts. 
Therefore, we find from direct FIR detections that the ULIRG 
fraction in LAEs decreases from z~2.4-2.8 to ~0.3, following 
the trend found for IR- detected galaxies between those redshifts 
(iLeFloc'het alll2005h . 

This possible evolution of the IR emission, along with the 



cates that either the Lya selection technique does not trace the 
same kind of objects with redshift, or the properties of Lya se- 
lected galaxies change with redshift. 

6.3. Dust attenuation 

The ratio between IR and UV luminosities is the b est way to 
obtain the dust attenuation in galaxies. Adopting the iBuat et al.l 
(2005) calibration, the dust attenuation for our PACS-detected 
Lya emitting galaxies can be obtained with the expression 



= -0.0495X 5 + 0.4718.x 2 + 0.8998.x + 0.2269 



(2) 



where, x = log {Lml L^uv)- The conversion from dust attenu- 
ation in NUV band into the co lor excess in t he ste llar contin- 
uum is carried out by using the lCalzetti et a l. (2000) reddening 
law. The values obtained are shown in Table Q] PACS-detected 
LAEs have high values of dust attenuation and, since they are 
selected via their IR emission, those values represent an up- 
per limit in dust attenuation of LAEs in the studied redshift 
range. This is a direct indication from FIR measurements that 
the Lya emission and the presence of dust are not mutually ex- 
clusive. This result could reinforce the scenario wher e the Lya 
emiss ion can be enhanced under a suitable geometry dNeufeldi 
1991), for which evidence has already been found at higher red- 
shifts dFinkelstein et all 120081). but not at o ur redshift or lower 
dGuaita et all 1201 U iFinkelstein et all 1201 lh . The dust attenua- 
tion derived from SED-fitting tends to be lower than obtained 
from direct IR/UV measurements. The ULIRG nature plays an 
important role in this underestimation, since for such as sources, 
the IR emission is more prominent than in other less IR luminous 
galaxies, and its inclusion in the analysis is essential to obtaining 
accurate results. 



6.4. SFR 

The combination of UV and IR luminosities provides the best 
determination of the SFR in galaxies. Assuming that all the 
light absorbed in the UV is in turn reradiated in the FIR, the 
total SFR of a galaxy, SFR, 0/fl ;, can be obtained as the sum of 
S¥Ruv,uncorrected plus a correctio n term associated to the dust 
emission in the FIR. Adopting the [Kennicuttl (1 1 998b calibration, 
that correction term can be written as 



SFR/stMoyr- 1 ] = 4.5 • 10- 44 L /R [ergs" 1 ] 



(3) 



where Ljr is defined in the same way than in Section loTTI In this 
way, SFRtotal can be inferred from 



SFR to , fl/ = SFR 



UV, uncorrected 



+ SFR„ . 



(4) 



reported evolution in the UV luminosi ty (Section 



Lya luminosity function from z~ 1 .0 dCowie et al 



M and in the 



201 lh . indi- 



The values of SFR, 0M / obtained for the PACS-detected galax- 
ies are also shown in TableQ] Such sources have high SFR, oto /, all 
above 2OOM yr~ 1 . PACS-deteced sources are the reddest ones of 
the sample and the fastest star formers. Although the number of 
such kind of galaxies is low, their SFRs can therefore be con- 
sidered as the maximum value of SFRs for LAEs in this redshift 
range. It can be also seen that the contribution of the SFR/« to 
the SFR, oto / is quite high, above 93% in all cases, and, there- 
fore, SED fitting with BC03 templates, which do not take the 
dust emission in the FIR into account, systematically underesti- 
mates this quantity. Since we do not have a statistically signifi- 
cant sample of LAEs detected in the FIR, we cannot evaluate a 
typical factor of underestimation for LAEs in our redshift range. 
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Table 2. Total infrared luminosities for the MIPS-24yum detected 
LAEs. PACS-detected LAEs, which are also detected in MIPS- 
24yum are not included in this table. 



Object 


Redshift 


Log(L //; /L Q ) 


GOODS_LRb_002_ 1 _q3_8_ 1 


2.296 


12.34 


GOODS.LRb_002_l_q2_60_l 


2.181 


11.30 


GOODSiRb_dec06_2.q2.12_l 


2.318 


12.10 


GOODS_LRb_dec06_3_q 1 _2 1 _ 1 


2.230 


11.97 


GOODS_LRb_002_l_q4_66_l 


2.586 


12.29 


GOODS_LRb_dec06_3_q4_53_l 


2.702 


11.62 



In lOteo et all (1201 lbh we found that, for LAEs at z~0.3, the con- 
tribution of the FIR emission is more than 60% in most cases, so 
that SFR based on SED fitting is underestimating the SFR by a 
factor greater than 2. According to the IR evolution reported in 
Section l6\2l it could be expected that the FIR contribution can be 
higher than 60% for a high percentage of the sample. Actually, 
LAEs at z~2.0 with L/r>10 105 L g are unlikely to be detected 
with the current FIR surveys, and their IR contribution to the 
SFR MM /, assuming a SFRyy.^^c^-lOMoyr 1 , is about 50%. 

By using the dust attenuation derived from SED fitting, we 
can correct the S¥Kuv, U ncorrected of the PACS-detected sources 
to estimate the SFR M/II /. Their values are also shown in Table 
[1] It can be seen that the SFR, OM / obtained in such a way are 
much lower than the values ob tained from the IR emission. This 
agrees wifh lWuvts et alj d201 lb . who find that SED fitting tend 
to underestimate the total amount of star formation in IR-bright 
galaxies at z>2.5. 

6.5. MIPS-24/jm counterparts of LAEs 

In addition to the analyzed measurements in PACS-160/vm, six 
other PACS-undetected galaxies of the total sample of 144 have 
MIPS-24yum detections in the FIR regime, which were found fol- 
lowing the same criterion as used for PACS counterparts. All 
these galaxies are classified as LAEs. None of them are within 
the MUSIC footprint, so we cannot analyze the distinctive prop- 
erties of MIPS-24/vm-detected LAEs regarding the rest-frame 
UV to mid-IR SED fitting. The MIPS-24jum detections also in- 
dicate the presence of IR-bright, i.e. dusty, LAEs. At 2.0<z<3.5, 
the MIPS-24//m band has a significative contribution of PAH 
molecule emission, which does not allow us to determi ne the Ljr 
for MIPS -detected LAEs. Actually, lElbaz et all d2010h study the 
validity of MIPS-24/mi extrapolations to the L/«, finding that 
they are valid up to z~1.5 and for objects that are below the 
ULIRG limit. At z> 1 .5, MIPS-24jum extrapolations tend to over- 
estimate the Ljr in a factor of about 2 or 3. Despite this overes- 
timation, we show the L/« luminosities for the MIPS-24/im de- 
tected LAEs in Table|2]as inferred from MIPS-24jum single band 
extrapolation, in a similar way to Sect. 16. II for PACS-detected 
sources. Owing to the PACS catalogs used in Section [6] PACS- 
detected LAEs are also detected in MIPS-24yum but we do not 
consider those objects here because they have been more ac- 
curately studied in previous sections. It can be seen in Table [2] 
that, although overestimated, MIPS-24/mi detected LAEs have 
high values of L/#, placing them in, at least, the LIRG regime, 
L/fl>10 11 L G . Therefore, in total we gather a sample of ten red 
and dusty LAEs, which are opposite to the classical idea of LAEs 
as dust-free galaxies. 



7. Summary and conclusions 

In this work we have analyzed the rest-frame UV to IR SED of a 
sample of 56 Lya emitting galaxies at 2.0<z<3.5. According to 
their Lya EW„ I( -/ ra/H we distinguish between LAEs and non- 
LAEs as those Lya emitting galaxies with Lya EW rej/ _/ rame 
above and below 20A, respectively. This value is the typical 
threshold in narrow-band searches. Our main conclusions are as 
follows. 

1. Individual SED fittings for the studied LAEs indicate that 
they are mostly young galaxies with ages below 100 Myr, 
although some of them have SEDs compatible with pop- 
ulations older than 300 Myr. The dust attenuation ranges 
from low, E(B-V)~0.0, to high values, E(B-V)~0.3, with a 
median of E(B-V) of 0.15. The SFRs, as derived from UV 
dust-uncorrected measurements, range from about 2 to 100 
M yr _1 , mostly having SFRf/v.^orrec/^-lOMoyr 1 . 

2. We find that LAEs at 2.0<z<3.5 show a wide range of prop- 
erties, rather than having a double nature as suggested in pre- 
vious works. This result has been possible to achieve because 
our LAEs are continuum-bright objects and, therefore, an in- 
dividual SED fitting for each source could be carried out. 

3. The variety of the physical properties of LAEs is also re- 
flected in their morphology. We find many different struc- 
tures, from bulge-like galaxies to highly clumpy systems or 
chain galaxies. Fitting the bulge-like LAEs to Sersic profiles, 
we found sizes distributed within 2.6 and 3.1 kpc, compati- 
ble with those reported in other works at higher, similar and 
lower redshifts. 

4. We find four Lya emitting galaxies with PACS-FIR coun- 
terparts, two LAEs, and two non-LAEs. This indicates that 
some LAEs have such high dust content that their emis- 
sion can be directly detected. Their rest-frame UV to mid-IR 
SEDs are compatible with objects with high dust attenua- 
tion. FIR-detected LAEs are objects with high SFR, above 
2OOM yr~ 1 , and their L !R classify them as ULIRGs. LAEs 
with an ULIRG nature are not seen at z~0.3, suggesting that 
there is an evolution in their IR emission with redshift. This, 
along with the rapid evolution of the Lya luminosity function 
from z~0.3 to ~1.0 reported in previous works and the evo- 
lution in the Lj/y from z~0.3 to z>2.0 found here indicates 
that either the physical properties in the UV, optical, and IR 
of objects that are selected via their Lya emission change 
with redshift or else the Lya selection technique would seg- 
regate different types of galaxies at different redshifts. 

5. We obtain the dust attenuation in the PACS-detected LAEs 
with MUSIC counterparts by combining UV and FIR infor- 
mation, without the uncertainties of the SED-fitting based 
methods. The results obtained in this way show that they 
have color excesses greater than 0.30 in their stellar con- 
tinuum. This indicates that PACS-detected LAEs are dusty 
objects, despite having a Lya line emission in their spectra. 
This confirms from direct FIR observations that dust and Lya 
emission are not mutually exclusive. 
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Table 1. Dust attenuation as derived from SED fitting, E(B-V). S [BC03], and from the IR/UV ratio, E(B-V) S [IR/UV], dust- 
uncorrected SFR, SFR uv , U ncorr, dust-corrected SFR, SFRd ust - corr , and total SFR, SFR foto ;, for the four PACS-detected sources. 



Object 


E(B-V)., [BC03] 


E(B-V)., [IR/UV] 


SFR uv ,imcorr [M Q Vr '] 


SFR dust - corr [M yr '] 


SFR toto/ [M yr-'] 


GOODS_LRb_00 1 _q 1 _8_ 1 


0.20 


0.40 


19.6 


74.9 


242.3 


GOODS_LRb_dec06_3_q3_60_l 


0.4 


0.48 


11.9 


172.7 


311.2 


GOODS_LRb_00 1 _q3_9_ 1 


0.30 


0.57 


10.7 


80.1 


452.8 


GOODS_LRb_001_q2_9_l 


0.15 


0.42 


25.6 


70.3 


341.5 



a contract to ESA managed by the Herschel/Planck Project team by an indus- 
trial consortium under the overall responsibility of the prime contractor Thales 
Alenia Space (Cannes), and including Astrium (Friedrichshafen) responsible for 
the payload module and for system testing at spacecraft level, Thales Alenia 
Space (Turin) responsible for the service module, and Astrium (Toulouse) re- 
sponsible for the telescope, with in excess of a hundred subcontractors. 



References 

Balestra, I., Mainieri, V., Popesso, P., et al. 2010, A&A, 512, A12+ 
Bond, N., Gawiser, E., Guaita, L., et al. 201 1, ArXiv e-prints 
Bond, N. A., Gawiser, E., Gronwall, C, et al. 2009, ApJ, 705, 639 
Bongiovanni, E., Oteo, I., Cepa, J., et al. 2010, A&A, 519, L4+ 
Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 
Buat, V., Iglesias-Paramo, J., Seibert, M., et al. 2005, ApJ, 619, L51 
Calzetti, D., Armus, L., Bohlin, R. C, et al. 2000, ApJ, 533, 682 
Chapman, S. C, Blain, A. W., Smail, I., & Ivison, R. J. 2005, ApJ, 622, 772 
Chary, R. & Elbaz, D. 2001, ApJ, 556, 562 
Cowie, L. L., Barger, A. J., & Hu, E. M. 2010, ApJ, 711, 928 
Cowie, L. L., Barger, A. J., & Hu, E. M. 2011, ApJ, 738, 136 
Cowie, L. L. & Hu, E. M. 1998, AJ, 115, 1319 
Deharveng, J., Small, T., Barlow, T. A., et al. 2008, ApJ, 680, 1072 
Elbaz, D., Dickinson, M., Hwang, H. S., et al. 2011, A&A, 533, A119+ 
Elbaz, D., Hwang, H. S., Magnelli, B., et al. 2010, A&A, 518, L29+ 
Elmegreen, B. G., Elmegreen, D. M., Fernandez, M. X., & Lemonias, J. J. 2009, 
ApJ, 692, 12 

Feldmann, R., Carollo, C. M, Porciani, C, et al. 2006, MNRAS, 372, 565 
Finkelstein, S. L., Cohen, S. H, Malhotra, S., & Rhoads, J. E. 2009a, ApJ, 700, 
276 

Finkelstein, S. L., Cohen, S. H, Moustakas, J., et al. 2011, ApJ, 733, 117 
Finkelstein, S. L., Rhoads, J. E., Malhotra, S., & Grogin, N. 2009b, ApJ, 691, 
465 

Finkelstein, S. L., Rhoads, J. E., Malhotra, S., Grogin, N., & Wang, J. 2008, ApJ, 
678,655 

Gawiser, E., Francke, H, Lai, K., et al. 2007, ApJ, 671, 278 

Gawiser, E., van Dokkum, P. G., Gronwall, C, et al. 2006, ApJ, 642, L13 

Grazian, A., Fontana, A., de Santis, C, et al. 2006, A&A, 449, 951 

Gronwall, C, Ciardullo, R., Hickey, T., et al. 2007, ApJ, 667, 79 

Guaita, L., Acquaviva, V., Padilla, N., et al. 2011, ApJ, 733, 114 

Guaita, L., Gawiser, E., Padilla, N., et al. 2010, ApJ, 714, 255 

Hogg, D. W., Cohen, J. G., Blandford, R., & Pahre, M. A. 1998, ApJ, 504, 622 

Hopkins, A. M. 2004, ApJ, 615, 209 

Hopkins, A. M. & Beacom, J. F. 2006, ApJ, 651, 142 

Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189 

Kornei, K. A., Shapley, A. E., Erb, D. K., et al. 2010, ApJ, 711, 693 
Lai, K., Huang, J., Fazio, G., et al. 2008, ApJ, 674, 70 
Le Floc'h, E., Papovich, C, Dole, H, et al. 2005, ApJ, 632, 169 
Lutz, D., Poglitsch, A., Altieri, B., et al. 2011, A&A, 532, A90+ 
Madau, P. 1995, ApJ, 441, 18 

Murayama, T., Taniguchi, Y., Scoville, N. Z., et al. 2007, ApJS, 172, 523 

Neufeld, D. A. 1991, ApJ, 370, L85 

Nilsson, K. K. & M0ller, P. 2009, A&A, 508, L21 

Nilsson, K. K., M0ller, P., MSller, O., et al. 2007, A&A, 471, 71 

Nilsson, K. K., Ostlin, G., M0ller, P., et al. 2011, A&A, 529, A9+ 

Nilsson, K. K., Tapken, C, M0ller, P., et al. 2009, A&A, 498, 13 

Nordon, R., Lutz, D., Shao, L., et al. 2010, A&A, 518, L24+ 

Oke, J. B. & Gunn, J. E. 1983, ApJ, 266, 713 

Ono, Y., Ouchi, M., Shimasaku, K., et al. 2010, MNRAS, 402, 1580 
Oteo, I., Bongiovanni, A., Perez Garcia, A. M., et al. 2011a, ApJ, 735, L15+ 
Oteo, I., Bongiovanni, A., Perez Garcia, A. M., et al. 2011b, Submitted 
Ouchi, M., Shimasaku, K., Akiyama, M., et al. 2008, ApJS, 176, 301 
Ouchi, M., Shimasaku, K., Furusawa, H, et al. 2010, ApJ, 723, 869 
Peng, C. Y, Ho, L. C, Impey, C. D., & Rix, H. 2010, AJ, 139, 2097 
Pentericci, L., Grazian, A., Fontana, A., et al. 2007, A&A, 471, 433 
Pentericci, L., Grazian, A., Scarlata, C, et al. 2010, A&A, 514, A64+ 



Perez-Gonzalez, P. G., Rieke, G. H, Egami, E., et al. 2005, ApJ, 630, 82 
Pilbratt, G. L., Riedinger, J. R., Passvogel, T, et al. 2010, A&A, 518, L1 + 
Pirzkal, N, Malhotra, S., Rhoads, J. E., & Xu, C. 2007, ApJ, 667, 49 
Poglitsch, A., Waelkens, C, Geis, N., et al. 2010, A&A, 518, L2+ 
Popesso, P., Dickinson, M., Nonino, M., et al. 2009, A&A, 494, 443 
Reddy, N. A., Erb, D. K., Pettini, M., Steidel, C. C, & Shapley, A. E. 2010, ApJ, 
712, 1070 

Rhoads, J. E., Malhotra, S., Dey, A., et al. 2000, ApJ, 545, L85 
Salpeter, E. E. 1955, ApJ, 121, 161 

Santini, P., Fontana, A., Grazian, A., et al. 2009, A&A, 504, 751 
Sersic, J. L. 1968, Atlas de galaxias australes, ed. Sersic, J. L. 
Shioya, Y, Taniguchi, Y, Sasaki, S. S., et al. 2009, ApJ, 696, 546 
Wuyts, S., Forster Schreiber, N. M., Lutz, D., et al. 2011, ApJ, 738, 106 



9 



I. Oteo et al.: Stellar populations and Herschel FIR-counterparts of hya emitters at 2 < z < 3.5 in GOODS-South 





IxlO 4 




GOODS_LRb_001_ql_8_l 






IxlO 3 










IxlO 2 

IxlO 1 

1x10° 

IxlO" 1 
1 vi n~ 2 

1 X 1 u 






- 






' 1 1 1— 


- 
- 




IxlO 4 




GOODS_LRb_dec06_3_q3_60_l 






IxlO 3 










IxlO 2 










IxlO 1 
1x10° 
IxlO -1 

1 v 1 A -2 

1 X 1 U 




' 1 1 1— <-w 


- 

- 




IxlO 4 




nOOnS T Rh D01 Q 1 


1— 




IxlO 3 








i 1 


IxlO 2 










IxlO 1 




of 






1x10° 






- 




IxlO" 1 






- 




1 vlfl" 2 

1 X 1 u 




1 1 1 1— w 






IxlO 4 




GOODS_LRb_001_q2_9_l 


1— 




IxlO 3 










IxlO 2 






\ — 




IxlO 1 




_ « r - T -CX56"cM-' / 






1x10° 




^-oco-o^o 






lxicr 1 










IxlO" 2 









10 4 10 5 10 6 10 7 

Wavelength [Angstroms] 



Fig. 6. Optical-to-FIR observed SED of the PACS-detected 
galaxies. Blue dots represent the multi-wavelength photometry 
from U-band to MIPS-24-yum, and red dots are the PACS-100/mi 
and/or PACS-160yum fluxes. Black curve represents the combi- 
nation of BC03 and CE01 templates that best fit the photometric 
data of each galaxy. 
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Appendix A: Physical properties of the studied Lya 
emitting galaxies 
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Table A.l. Physical properties of our sample of 56 galaxies at 2.0<z<3.5 as derived from SED fitting with BC03 templates. 



Name 


Redshift 


Age [Myr] 


E(B-V)[BC03] 


SFK[jv,uncorrecled [M Q yr '] 


log(M»/M Q ) 


GOODS_LRb_dec06_3_q2_33_2 


2.0214 


29 ± 34 


0.25 ± 0.07 


6.7 


8.2 


GOODS_LRb_002_q3_86_l 


2.0526 


21+6 


0.30 ± 0.05 


17.0 


8.5 


GOODS_LRb_002_l_q3_88.1 


2.0600 


47 ± 8 


0.15 ±0.05 


4.8 


8.3 


GOODS_LRb_002_ 1 _q4_79_ 1 


2.1449 


25 ±6 


0.15 ±0.23 


20.0 


8.6 


GOODS _LRb_00 1 _q3_24.2 


2.1659 


15 + 6 


0.20 ± 0.06 


9.1 


8.1 


GOODS_LRb_dec06_ 1 .q2_4 1 _1 


2.1809 


541 ±6 


0.10 ± 0.05 


2.3 


9.1 


GOODS_LRb_dec06_l.q4_15_2 


2.2128 


121 ± 17 


0.10 ±0.05 


20.4 


9.3 


GOODS_LRb_001_q2_7_l 


2.2158 


383 + 10 


0.15 ±0.05 


2.5 


8.9 


GOODS_LRb_dec06_l_q2_48_l 


2.2991 


11 ±32 


0.20 ± 0.09 


3.2 


7.5 


GOODS _LRb_002_ l_ql_16_2 


2.3115 


5±5 


0.25 ±0.19 


15.9 


7.9 


GOODS_LRb_002_l_q4_29_l 


2.3143 


11 ±24 


0.30 ± 0.06 


9.7 


8.0 


GOODS_LRb_dec06_3_q3_5 1_1 


2.3166 


305 ±9 


0.15 ±0.05 


6.2 


9.2 


GOODS_LRb_dec06_3_q2_22_l 


2.3170 


31 ±21 


0.15 ±0.05 


17.3 


8.7 


GOODS_LRb_dec06_ l_ql_ll_3 


2.3200 


55 ±22 


0.20 ± 0.06 


27.5 


9.1 


GOODSiRb_dec06_l.ql.15_l 


2.4000 


121 ± 5 


0.00 ± 0.05 


5.4 


8.8 


GOODS_LRb_001_ql_8_l 


2.4284 


1750 ± 13 


0.20 ±0.14 


19.6 


10.5 


GOODS_LRb_002_l_q4_69_l 


2.4304 


33 ±32 


0.25 ± 0.05 


13.1 


8.6 


GOODS_LRb_dec06_3_q3_53_3 


2.4645 


39 ±5 


0.10 ± 0.05 


4.5 


8.2 


GOODS_LRb_dec06_l_q3_60_l 


2.4834 


13 ±6 


0.15 ±0.05 


10.9 


8.1 


GOODS_LRb_002_l_q4_l_l 


2.4835 


121 ±7 


0.05 ± 0.05 


22.1 


9.4 


GOODS_LRb_dec06_3_q3_60_l 


2.5144 


9 ± 18 


0.40 ± 0.05 


11.9 


8.0 


GOODS _LRb_002_ 1 _q 1 _26_ 1 


2.5600 


1 ± 5 


0.40 ±0.16 


13.4 


7.1 


GOODS_LRb_00 1 _q3_47_2 


2.5641 


383 ±9 


0.00 ± 0.05 


4.8 


9.2 


GOODS_LRb_001_q2_14_l 


2.5671 


215 ±9 


0.15 ±0.05 


9.2 


9.2 


GOODS_LRb_002_l_q4_32_l 


2.6159 


305 ± 8 


0.00 ± 0.05 


7.6 


9.3 


GOODS_LRb_002_l_ql_62_l 


2.6191 


341 ± 8 


0.00 ± 0.05 


7.8 


9.4 


GOODS_LRb_001_q2_37_l 


2.6489 


77 ± 10 


0.10 ±0.05 


12.5 


8.9 


GOODS_LRb_dec06_l_q3_39_2 


2.6607 


5±6 


0.25 ± 0.07 


11.4 


7.7 


GOODS_LRb_dec06_l_q3_32_l 


2.6692 


153 ± 13 


0.10 ± 0.05 


8.0 


9.0 


GOODS_LRb_00 1 _q3_69_ 1 


2.6796 


7 ±4 


0.10 ±0.05 


4.6 


7.5 


GOODS_LRb_00 1 _q2_3 1 _ 1 


2.6914 


1210 ± 10 


0.10 ±0.05 


9.1 


10.0 


GOODS_LRb_00 1 _q3_60_ 1 


2.6933 


171 ± 11 


0.20 ± 0.05 


7.2 


9.0 


GOODS_LRb_dec06_2_q3_44_l 


2.7034 


305 ± 14 


0.15 ±0.05 


16.0 


9.6 


GOODS_LRb_00 1 _q2_39_ 1 


2.7184 


101 ±6 


0.05 ± 0.05 


6.4 


8.8 


GOODS_LRb_00 1 _q2_44_2 


2.7237 


3±2 


0.30 ±0.19 


104.6 


8.4 


GOODS_LRb_dec06_l_q4_5_2 


2.7269 


3±5 


0.30 ± 0.08 


10.7 


7.5 


GOODS_LRb_001_q3_9_l 


2.8079 


27 ± 10 


0.30 ± 0.05 


10.7 


8.4 


GOODS_LRb_001_q2_9_l 


2.8121 


61 ±9 


0.15 ±0.05 


25.6 


9.1 


GOODS_LRb_dec06_2_q3_5 1_1 


2.8130 


101 ± 10 


0.10 ± 0.05 


9.1 


8.9 


GOODS_LRb_001_q4_22_l 


2.9659 


429 ±3 


0.00 ± 0.05 


3.2 


9.1 


GOODS_LRb_001_q2_13_l 


3.0058 


5±2 


0.25 ±0.19 


119.2 


8.7 


GOODS_LRb_00 1 _q3_48_ 1 


3.0201 


9±6 


0.30 ± 0.10 


16.4 


8.1 


GOODS_LRb_00 1 _q3_89_2 


3.0300 


121 ± 11 


0.15 ±0.05 


10.6 


9.1 


GOODS_LRb_dec06_l_q4_7_l 


3.0872 


21 ±27 


0.20 ± 0.07 


12.3 


8.4 


GOODS_LRb_002_l_ql_52_2 


3.0920 


39 ±8 


0.10 ±0.05 


3.9 


8.1 


GOODS_LRb_dec06_l_q4_45_3 


3.1316 


121 ±9 


0.10 ± 0.05 


6.7 


8.9 


GOODS_LRb_001_q2_84_l 


3.1677 


9±6 


0.15 ±0.05 


12.1 


8.0 


GOODS_LRb_00 1 _q3_84_ 1 


3.1733 


763 ± 10 


0.05 ± 0.05 


6.4 


9.6 


GOODS_LRb_dec06_3_q3_42_l 


3.2000 


193 ± 18 


0.10 ±0.05 


6.8 


9.1 


GOODS_LRb_dec06_ 1 _q3_4 1 _1 


3.2078 


69 ± 10 


0.10 ±0.05 


27.0 


9.2 


GOODS_LRb_dec06_3_q3_65_2 


3.2293 


57 ± 16 


0.20 ± 0.05 


21.5 


9.0 


GOODS_LRb_00 1 _q3_7 1 _2 


3.2479 


87 ±25 


0.15 ±0.06 


9.3 


8.9 


GOODS_LRb_002_l_ql_33_l 


3.3565 


5 ±4 


0.20 ± 0.05 


8.7 


7.6 


GOODS_LRb_001_q2_30_l 


3.3850 


153 ±6 


0.05 ± 0.05 


8.1 


9.0 


GOODS_LRb_002_l_q2_71_l 


3.4724 


25 ±9 


0.15 ±0.05 


21.5 


8.7 


GOODS_LRb_dec06_3_q3_3 1_1 


3.5400 


37 ±6 


0.10 ±0.05 


4.6 


8.2 



